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Abstract. Without amplification, magnetic fields in expanding ejecta of young supernova remnants (SNRs) will be 
orders of magnitude below those required to shock accelerate thermal electrons, or ions, to reiativistic energies or to 
produce radio synchrotron emission at the reverse shock. The reported observations of such emission give support 
to the idea that diffusive shock acceleration (DSA) can amplify magnetic fields by large factors. Furthermore, the 
uncertain character of the amplification process leaves open the possibility that ejecta fields, while large enough 
to support radio emission and DSA, may be much lower than typical interstellar medium values. We show that 
DSA in such low reverse shock fields is extremely nonlinear and efficient in the production of cosmic-ray (CR) 
ions, although CRs greatly in excess of mc 2 are not produced. These nonlinear effects, which occur at the forward 
shock as well, are manifested most importantly in shock compression ratios ^> 4 and cause the interaction region 
between the forward and reverse shocks to become narrower, denser, and cooler than would be the case if efficient 
cosmic-ray production did not occur. The changes in the SNR structure and evolution should be clearly observable, 
if present, and they convey important information on the nature of DSA and magnetic field amplification with 
broad astrophysical implications. 
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1. Introduction 

It is clear in many supernova remnants (SNRs) that the 
forward, blast wave shock, interacting with the interstellar 
medium (ISM) magnetic field, produces radio (and some- 
times X-ray) synchrotron emission. Presumably this is ac- 
complished when the forward shock accelerates thermal 
and pre-existing cosmic-ray (CR) electrons by diffusive 
shock acceleration (DSA) (also called the first-order Fermi 
mechanism). The reverse shock, however, will not produce 
reiativistic electrons from thermal ones if the only mag- 
netic field that is present is the ambient field from the pro- 
genitor star. Any progenitor field will be vastly diluted by 
expansion and flux freezing and, for expected white dwarf 
or massive star magnetic field values, after < 100 yr, will 
fall below levels necessary to support particle acceleration 
to radio emitting energies. 

For example, if the surface field of a white dwarf of 
radius Rwd — 10 7 m is i%D = 10 9 G, the diluted mag- 
netic field after ~ 100 yr at the reverse shock, which is 
typically at a radius i? s k ~ 1 pc from the explosion site, 
is B ~ B WD {R WD /R sk ) 2 ~ 10~ 10 G. For lO" 10 G, the 
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diffusive acceleration time to 10 GeV is <; 100 yr and the 
upstream diffusion length at a 10 4 km s _1 shock of a 10 
GeV electron is > 1 pc, making the production of radio 
emitting electrons unlikely. Similar results are expected 
for massive progenitors. Furthermore, with such low fields, 
even if reiativistic electrons are present in large numbers 
from compressed pre-existing cosmic rays or whatever, the 
radio emissivity may be too low to be observable. 

The expanded ejecta bubble may be one of the lowest 
magnetic field regions in existence and if reverse shocks 
in some SNRs are accelerating electrons by DSA to ra- 
dio emitting energies or higher (as has bee n suggested by 
recent radio and X-ray observations, e.g., iGotthelf et all 
l200lUDeLanev et al.112002. iRho etHhoO^ . there are im- 
portant consequences for: 

— magnetic field generation and amplification in strong 
shocks, 

— the intrinsic efficiency of DSA and cosmic-ray produc- 
tion, including heavy elements, and 

— the structure and evolution of the X-ray emitting in- 
teraction region between the forward (FS) and reverse 
shocks (RS). 
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At the fo rward shock in some SNRs (see, in partic- 
ular, Cas A; IVink fc Lamind 12003). there is convincing 
evidence for magnetic fields far greater than normal 
ISM values. It has been suggested that the diffusive 
shock acceleration proces s can amplify ambie nt fields 
(e.g., iLucek fc BelH l2000t iBell fc Lucekl l200ll) to the 
observed levels, an d evidence for this is mounting in spe- 
cific SNRs, fe.g.. iBerezhko Ksenofontov. fc VolM 120031: 



Berezhko. Piihlhofer. fc Void 1200.4 iBerezhko fc Volkl 
If such amplification occurs generally, SNRs will 
be ca pable of accelerating co s mic r a ys to above 10 17 eV 
(e.g., IPtuskin fc Zirakashvilil 120031: iDrurv et all l2003|) . 
possibly solving the decades old problem of smoothly 
generating CRs to the spectral knee near 10 15 eV and 
beyond. Since DSA is expected to occur in diverse 
environments on all astrophysical scales, the confirmation 
and characterization of magnetic field amplification is of 
extreme importance. 

In this paper, we first discuss in §2 the general effects 
of a low magnetic field on the efficiency of DSA. We show 
that nonlinear effects, most noticeable in producing com- 
pression ratios 3> 4, can be extremely large for a range 
of magnetic field below that of the average ISM (i.e., for 
B < Sism - 3 x 10~ 6 G). We then discuss in §3 the rel- 
evance of low fields and strong nonlinear effects in DSA 
at reverse shocks in SNRs. To our knowledge, this is the 
first attempt to consider such effects in an evolutionary 
model of SNRs. A number of aspects concerning the ac- 
celeration process, the nature of the magnetic field, and 
the physical conditions in the unshocked ejecta material 
(e.g., temperature, ionization fraction, etc.) are not well- 
known. Because of these uncertainties, we show a number 
of examples where important parameters are varied over 
fairly wide ranges. We emphasize, however, that the ef- 
fects of efficient DSA on the structure and evolution of 
SNRs may be large and current instruments should be 
sensitive enough to importantly constrain many of these 
poorly known parameters. 

2. Nonlinear Diffusive Shock Acceleration 

Nonlinear DSA is a complex process that is difficult to 
describe completely. In order to allow the coupling of 
the particle acceleration to a hydrodynamic model of 
SNR evolution, we use an approximate, alg e braic m odel 
of DSA dev e loped by IBerezhko fc Ellisonl l|l999l) and 
lEllison et alJ ll2000l) Whil e more complete models exist 
fe.g.. IBerezhko et alJ ll99fi). they tend to be more compu- 
tationally intensive and not as easy to include in a hydro- 
dynamical simulation. Despite the simplifications made in 
our acceleration model, we believe it adequately describes 
the essential physics when the maximum momentum of ac- 
celerated particles p max ^> Am p c. It is less accurate when 
Pmax ~ ArripC as we describe more fully below. Here, A 
is the mass number, m p is the proton mass, and c is the 
speed of light. 

In a complicated, nonlinear fashion, the acceleration 
efficiency (i.e., the fraction of total ram kinetic energy go- 



ing into superthermal particles) 1 depends on the sonic 
and Alfven Mach numbers (Ms = \J PoV^/ilPo) and 
Ma — \/^PoV^/B respectively), on the particle injec- 
tion and on the maximum particle momentum achieved 
Pmax- Here, po is the unshocked mass density, Pq is the 
unshocked pressure, B is the unshocked magnetic field, 
T4k is the shock speed, and 7 is the ratio of specific 
heats. The particle injection is modeled by two param- 
eters: the injection efficiency 7yi„j, i.e., the fraction of total 
particles which end up with superthermal energies, and 
Ai n j, which determines the value of the injection momen- 
tum (= Ai n jj4m p C S 2, where C S 2 is the sound speed in 
the downstream region). The parameter Ai n j, by defini- 
tion >JjJs_ai^Wtranlxtaken to be 4 in our calculations 
(see IBerezhko fc Ellisonl! 19991 for a full discussion). 

For simplicity, in all our examples we assume a single 
ion species, generally protons, but in §4.4, oxygen, with 
the electron temperature equal to the ion temperature 
= 10 4 K. We ignore any wave damping from neutral ma- 
terial (see §5.2.5 below for a discussion). The maximum 
momentum protons achieve is determined by setting the 
acceleration time equal to the SNR age isNR, or by setting 
the diffusion length of the highest energy particles equal to 
some fraction, / s k, of the shock radi us i? s i< , whichever g ives 
the lowest p max (see, for example, iBaring et al.lll999|) . 

We assume strong turbulence (i.e., Bohm diffusion) so 
that the scattering mean free path is on the order of the 
gyroradius, i.e., A ~ rj m { p r g , with ?7 m fp = 1. The magnetic 
field strength is thus an important factor in determining 
Pmax- If the turbulence is, in fact, weaker than Bohm dif- 
fusion (r) m f p > 1), Pmax will be smaller for a given back- 
ground B and the shock compression ratio rtot will be less. 
Consistent with assuming strong turbulence (?7 m fp = 1), 
we set, following I Volk et, alJ l|2002h . the downstream mag- 
netic field £?2 



2rl t /3 Bo 2 



In the acceleration model of IBerezhko fc Ellisonl 
the magnetic field also enters in a calculation of 
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the transfer of energy from energetic particles to the back- 
ground gas via Alfven waves, i.e., 



up ' 
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CR 
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where Va = B/y/Anp is the Alfven velocity, u is the flow 
speed, Pg a s is the pressure in the background gas, and 
Pgr is the pressure in the relativistic particles, i.e., cos- 
mic rays. It is implicitly assumed that the turbulence satu- 
rates when SB /B ~ 1 and the wave energy is then rapidly 
damped to heat. Thus, magnetic field amplification is not 
included in this description. 

It is via this energy transfer from energetic particles 
that the magnetic field has its largest effect on the accel- 
eration process and even small amounts of background 



1 Note the difference between injection efficiency 77inj and 
acceleration efficiency. The acceleration efficiency can and does 
vary with shock parameters even for a constant ryinj. 

2 Everywhere, the subscript (2) implies unshocked 
(shocked) quantities. 
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heating from the damping of Alfven waves can signifi- 
cantly reduce the acceleration efficiency compared to the 
case where only adiabatic heating is included. 

As explained in detail in iBerezhko fe Ellison! (^999), 
compression ratios > 4 occur in DSA for two reasons. 
First, as relativistic particles are produced and contribute 
significantly to the total pressure, their softer equation 
of state makes the shocked plasma more compressible (as 
7 — > 4/3, r to t — > 7). Second, as the highest energy particles 
escape from the shock, they drain away energy flux which 
must be compensated for by ramping up the overall com- 
pression ratio to conserve the fluxes. Just as in radiative 
shocks, this is equivalent to 7 — > 1 and r to t can become ar- 
bitrarily large. Compression ratios as large as we show here 
occur in independent steady-state calculati ons of nonlin- 
ear DSA which account for particle loss fe.g..lEichlerll984l 
I.Tones fc Ellisoi3ll99ll lMalkovlll99i lBlasill200^ . We note 
that as the overall compression ratio increases (r to t > 4), 
the subshock compression ratio, r su b, which is responsi- 
ble for heating the gas, becomes less than the test-particle 
value (r su b < 4), causing the temperature of the shocked 
gas to drop below test-particle values. 

Two important qualifications must be made concern- 
ing particle escape and the production of compression ra- 
tios greater than 7. First, the explicit assumption in the 
model of Berezhko and Ellison (and the others mentioned 
above) is that steady-state conditions apply and parti- 
cle acceleration is terminated as particles diffuse away 
from the shock. If, instead, the acceleration time be- 
comes comparable to the shock age while the diffusion 
length is still a small fraction of the shock radius, the 
acceleratio n process m ay terminate without particle es- 
cape fe.g.. lDrurvlll983h . It is less clear what happens in 
this case, although the work of Berezhko a nd co-workers 
(e.g.. lBerezhkolll996l IBerezhko et al.lll996j) suggests that 
for the forward shock in SNRs, geometrical factors deter- 
mine p m ax over most of the lifetime. These geometrical 
factors, i.e., the diffusion of particles upstream from the 
shock, the increase in the shock size and upstream vol- 
ume, the slowing of the shock speed, and the adiabatic 
cooling of the energetic particles, produce effects similar 
to those from escaping particles even though particles, in 
fact, remain in the SNR system (see a comparison of the 
modeling of SN1006 using the simple, steady-state model 
used here wit h the time-depende nt model of Berezhko and 
co-workers in lEllison et ai1l2000j) . 

The situation is more uncertain for the inward fac- 
ing, reverse shock since the upstream region has a finite 
volume and particles streaming far upstream can conceiv- 
ably reach the shock on the opposite side of the explosion 
site without being lost. In this case, particles may still be 
lost from spatial effects if they diffuse far enough down- 
stream to reach the contact discontinuity. However, large 
magnetic fields are expected at the contact discontinuity 
due to compression and stretching of the field lines by 
Ray leigh- Taylor instabilities, and might act as a magnetic 
wall confining the particles in the ejecta material. 



In any case, for either the forward or reverse shock, a 
precise determination of the compression ratio in a time- 
dependent situation requires a detailed knowledge of the 
wave-particle interactions in the self-generated magnetic 
turbulence responsible for particle diffusion. This knowl- 
edge does not yet exist so approximations remain neces- 
sary. We simply assume that particles at either the forward 
or reverse shock leave the system when their acceleration 
time > tsNR or their upstream diffusion length > / S k-R s k- 
In all of the SNR models we show in this paper, we arbi- 
trarily take f s k — 0.05 and note that p max and r to t would 
increase with increasing / s k- 

The second qualification is that the approximations 
in the Berezhko and Ellison model assume that p max ^> 
AnipC and thus, that a sizable fraction of the total pres- 
sure is in relativistic particles. Despite this limitation, 
we show cases where p max ^ 2Am p c and warn that our 
lower p max results have a greater intrinsic error than 
those with p max S> Am p c. In fact, from comparisons 
with Monte Carlo results (not shown here), we find that 
the Berezhko and Ellison model underestimates r to t when 

Pmax ~ AlTlpC. 

3. Efficient Diffusive Shock Acceleration in Weak 
Magnetic Fields 

In this section, we investigate how the acceleration effi- 
ciency depends on the upstream magnetic field Bq for a 
fixed set of parameters; i? s k, / s k, tsNRj and proton num- 
ber density n p . We chose a value of £snr large enough to 
ensure that p max is determined by diffusive escape. In this 
case, since we set r) m [ p = 1 and i? s k and fsk are fixed, p max 
depends only on Bq and V s k- 

Pmax — 3Ze/ sk i? sk £ Wc , (2) 

where Z is the charge number or ionization state and e is 
the electronic charge. 

Figure ^ shows the variation of the overall shock com- 
pression ratio rtot versus magnetic field (top panel) for 
different values of the injection parameter ryj n j with V^k 
set to 5000 km s^ 1 and (bottom panel) for different val- 
ues of the shock velocity with r/i n j set to 1CP 3 . In both 
panels, R s k = 3 pc, isNR = 1200 yr, n p — 0.3 cm -3 , and 
the plasma is composed of protons and electrons only. 

Each value of rj m ^ has two curves, one with only adi- 
abatic heating in the precursor and one with Alfven- 
wave heating in the precursor. The transfer of energy 
from the energetic particles to the background gas via 
Alfven waves reduces rtot dramatically as Bq increases. 
With Alfven-wave heating, very large compression ratios 
are obtained only for values of the magnetic field much 
lower than that typical of the interstellar medium (i.e., 
for Bo « 3 x 10~ 6 G). This compression ratio increases 
with decreasing magnetic field up to a maximum, whose 
position and intensity depend on the injection parameter, 
and then decreases to the test-particle value of 4. 

In the top panel, the sonic Mach number is fixed and 
the only relevant varying parameters are p m ax and the 



4 



Ellison et al.: Nonlinear Particle Acceleration at Reverse Shocks in Supernova Remnants 



150 



100 



50 



o 



c 

</) 

CD 



o 
o 







100 



50 - 







V sk =5000 km s' 



-1 - 




10 -9 -8 -7 -6 -5 -4 
Log 1Q B [G] 

Fig. 1. Compression ratio versus unshocked magnetic field 
for a constant shock speed V^k and different values of the 
injection efficiency 7/inj (top panel), and for a constant in- 
jection efficiency and different values of the shock veloc- 
ity (bottom panel). For each set of parameters, the upper 
curve corresponds to adiabatic heating in the precursor 
and the lower one to Alfven-wave heating in the precur- 
sor. In the latter case, r tot is dramatically reduced as B 
increases. 



Alfven Mach number, Ma- As Bo decreases, p max de- 
creases and Ma increases. For i?o<;2xlO~ 8 G, the mag- 
netic field is strong enough to allow a significant trans- 
fer of energetic particle energy through Alfven waves to 
heat, lowering the acceleration efficiency. The greater J5o, 
the smaller Ma and the more important this effect be- 
comes, causing r tot — * 4. For B < 2 x 10~ 8 G, Ma is 
large enough that it is no longer important, but now p max 
becomes small enough that the fraction of pressure in rel- 
ativistic particles drops below that required to maintain a 
strongly modified shock. As Bo decreases, the transition 
from a strongly modified shock with r to t 3> 4 to an un- 
modified one with rtot ~ 4 occurs and can b e extre mely 
abrupt. As explained in iBerezhko fc Ellisonl |l9991. the 



larger ?7i n j is, the more difficult it is to have a high-Mach 
number, unmodified solution. This is the reason that the 
maximum value of the compression ratio increases, and 
that the position of the maximum shifts towards lower 
Bq, as Tftnj increases. 

In the bottom panel, both Ms and Ma vary. With 
Alfven-wave heating, as B decreases, the curves for V4k = 
1000 and 2000 km s _1 end before the compression ratio 
drops to 4: in these cases, p max becomes low enough (^ 
2 m p c) to invalidate the approximations of the Berezhko 
and Ellison model. 

As Bo increases and Ma drops, the damping effects of 
the magnetic field increase and cause rtot to drop toward 
4, regardless of the sonic Mach number. For values of mag- 
netic field near 10~ 8 G in the Alfven-wave heating curves, 
the compression ratio peaks strongly as the shock veloc- 
ity and, therefore, Ms increase. As Bq decreases below 
~ 10~ 8 G, the lowering of p max , and subsequent reduc- 
tion of pressure in relativistic particles, causes the transi- 
tion to unmodified solutions regardless of Ms . The mag- 
netic field where r tot is a maximum shifts slightly toward 
lower Bq as Ms increases. Using the Berezhko and Ellison 
model, an approximate expression for the magnetic field 
strength, B*, at the transition point between unmodified 
and strongly modified solutions can be derived in two 
regimes (for p max < or > 100^4m p c): 



uj B* 



( Xl /VTo) vx 



■T2 



4 2 



(3) 



if p max > 100 Am p c , 



or 



xlVsk 

2 



X2 n 



(4) 



if p max < 100 Am p c , 



where 



^sub - 1 



(5) 



27- 



(Vtot/Vsub)" 



X2 = 



3/sk-Z'ei? s k 
Ai] mlp m p 



(6) 



At the transiti on, rtot, can be estimate d as 0.65 M|^ 4 
and r sub = 4 jBcrezliko fc Ellison! fl999() . The break at 
100 Am p c mirrors the break in the Berezhko and Ellison 
model between a three-component power law when p max > 
100 AnipC and a two-component power law when p max < 
100 Am p c. 

Figure ^ illustrates that even if extremely high com- 
pression ratios are theoretically possible, normal ISM 
magnetic field values, i.e., Bism *S 3x 10~ 6 G, with Alfven- 
wave heating, are sufficiently high to limit compression 
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ratios to r tot ^ 20, regardless of how efficient the injec- 
tion is or how high the sonic Mach number is. The energy 
converted to magnetic turbulence and heat lowers the sub- 
shock Mach number and the overall acceleration efficiency. 
Also, when Ma is low, the speed of the magnetic scatter- 
ing centers in the fluid (assumed to be Alfven waves) can 
become comparable to the shock speed, lowering the effec- 
tive difference (for acceleration) between the downstream 
and upstream flow speeds. 3 The low magnetic fields in 
expanding supernova ejecta offer a unique possibility of 
seeing non-radiative shocks with r to t <5 20. 



4. Particle acceleration at reverse shocks in SNRs 

4.1. Assumptions and parameters for the CR-hydro 
simulation 

We model the effects of cosmic-ray acceleration on the 
evolution of a SNR using a one-dimensional, cosmic- 
ray hydrodynamic (CR-hydro) simulation as described in 
lEllison et all (|2004|) . 

We initialize the CR-hydro simulation at some time to 
after the explosion with a power-law ejecta density dis- 
tribution, p c j oc r~" (of index n = 7), combined with 
a constant density plateau region at small radii. The 
plateau is required to keep the total ejecta mass finite. 
Beyond the ejecta, we assume a uniform ISM mass den- 
sity pisM = 5 x 10 -25 g cm -3 corresponding to a proton 
number density of 0.3 cm -3 . A constant density ISM is 
more appropriate for a Type la supernova, whereas a Type 
II supernova is likely to explode in a pre-SN stellar wind 
with a, p oc r~ 2 density structure. The presence of a stel- 
lar wind will affect the density and temperature structure 
of the shocked ejecta and this, in turn, will change the 
quantitative aspects of particle acceleration at the reverse 
shock. Qualitatively, however, the effects we describe for 
a uniform ISM will be present in Type II supernovae as 
well. 

The hydrodynamical simulation only recognizes the 
matter density, but the acceleration model depends on the 
ion species and is limited to a single species, as we discuss 
below. We assume that the ejecta speed varies linearly 
with radius from zero to some maximum speed V^ ax . In 
the simulation, the initial maximum radius of the ejecta is 
set by the maximum ejecta speed and to. Thus, the early 
stages of the simulation will depend on VJ° j ax and to . As 
long as the total kinetic energy and ejecta mass stay the 
same, however, the later evolution of the SNR is indepen- 
dent of both V^ ax and to- 

In Figure [3 we show how the early evolution depends 
on V^ ax for a particular set of SNR parameters. For all of 
our other examples, except those indicated in Fig. [21 we 
take V°j ax = 3 x 10 4 km s" 1 ~ 0.1c. As long as to is earlier 
than the time when p max becomes greater than 2 Am p c, 
as is the case in Fig. El our results are independent of to- 

3 Note that even though we consider cases where the speed 
of the scattering centers is high, we neglect second-order Fermi 
acceleration. 



O 



O 
O 



CO 



to 

> 



100 
50 

20 
10 
5 



2 = 



10" =- 



1000 



°^ 100 h 



I 10 r 

CL 



1 



1 1 — i — i i i 1 1 1 1 1 — i — i i i 1 1 1 r 

V mix= 4 >< 104 km /§* 



3x1 4 km/s 




H — I I I I III 



2x1 CT km/s 
H — I I I I M il 



RS speed 




H — I I I I III 1 — I I I I III h 



/ 



_l I I I I I 1_ 



1 10 100 

l SNR M 

Fig. 2. All curves are for the reverse shock and use the 
same set of parameters with a constant magnetic field 
B e j = 3 x 10~ 8 G. The only variable is the maximum 
speed of the ejecta V^ ax which is 4 x 10 4 km s _1 for the 
solid curves, 3 x 10 4 km s _1 for the dotted curves, and 
2 x 10 4 km s" 1 for the dashed curves. For comparison, 
we show the results obtained analytically, where V^ ax is 
not limited, for the same set of parameters using modified 
Chevalier solutions. While all curves converge after several 
decades, the early evolution depends on the different hy- 
drodynamic initial conditions. The simulation curves are 
truncated when p max < 2m p c. 



For comparison, we also show in Fig. (dot-dashed 
curve) the results obtained analytically for the same 
set of parameters using modified Chevalier solutions 
( Dcc ourchelle et al . 2000). Good agreement is reached af- 
ter a few tens of years and, in fact, all of our simulation 
results, after a few decades, are independent of our start- 
ing conditions and consistent with analytic solutions, as 
long as the self-similar conditions required for the analytic 
solutions are valid. 
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4.2. Constant ejecta magnetic field 

Given the general behavior of DSA in low magnetic fields, 
we begin our study of SNRs by assuming a constant ejecta 
magnetic field upstream from the reverse shock. 

For the CR-hydro model, we use the following param- 
eters: supernova explosion kinetic energy E SN = 10 51 erg, 
ejecta mass M C j = 1.4 M , and Bi SM = 3 x 1CP 6 G for the 
upstream ISM field (this is B Q for the forward shock). For 
the acceleration calculation, we take ?7i n j = 10~ 3 here and 
in all following models for both the forward and reverse 
shocks. 4 

The top panel of Fig.[2]shows the variation of the com- 
pression ratio with SNR age for different values of the 
constant upstream ejecta magnetic field B c) (Bq for the 
reverse shock). For B CJ > 3 x 1CP 6 G, the compression 
ratio is smaller than 20 for a SNR age larger than 100 yr. 
For B Ci < 3 x 10~ 7 G, however, the reverse shock compres- 
sion can be extremely high: r to t ~ 60 at isNR ~ 100 yr for 
B C j = 3 x 10~ 8 G (dashed curve). The analytical results for 
-B e j = 3 x 10~ 7 G are shown as long-dashed curves in the 
top two panels. As explained for Fig. the early evolu- 
tion depends on V^ J ax but good agreement is obtained for 
higher or, in any case, after ~ 10 yr. The change in 
the slope in the simulation curves at ~ 400 yr corresponds 
to the passage of the reverse shock from the power-law en- 
velope into the plateau. Such a transition is not included 
in the analytic result which is only valid while the reverse 
shock remains in the power-law profile. 

Lowering B e j causes the particle gyroradius and ac- 
celeration time to increase so in a shock of a given size 
and age, p max decreases, as seen in the middle panel of 
Fig. El Below some minimum value, B ej will be too weak 
to allow the acceleration of particles to radio emitting en- 
ergies. This opens up the possibility that a range of B Ci 
may exist greater than the minimum value needed to pro- 
duce observable radio emission but less than - 3 x 10" 6 G 
so that the full nonlinear effects of efficient Fermi acceler- 
ation of ions occurs. 5 

The efficiency of the DSA process can be extremely 
high. In the bottom panel of Fig. we show the instanta- 
neous acceleration efficiency E e fi, defined as the fraction 
of incoming energy flux (in the shock rest frame) put into 
relativistic particles. After a few decades, all of the models 
show E c ff <; 0.9. The curves in the bottom panel on the 
right show the fraction of the supernova explosion energy 
Es~n put into relativistic particles, i.e., Eqk/^SN' The re- 
verse shocks are able to put ~ 10% of E$n into CRs after 



4 In an actual SNR, of course, the injection efficiency might 
vary with time, vary over the shock surface, or be different at 
the forward and reverse shock s (as in our model of Kepler's 
SNR: IPecourchelle et al.ll2000T) . If the actual ?7i n j is less than 
10 -3 , the nonlinear effects we show will be less dramatic. 

5 Note that we do not show solutions in Fig. [3] for p ma x < 
2m p c, i.e., below the limit of validity of the Berezhko and 
Ellison model. 




Fig. 3. In each panel, RS results are shown for B C j = 3 x 
10~ 6 G (solid curves), 3 x 10~ 7 G (dot-dashed curves), 3x 
10~ 8 G (dashed curves), and 3 x 10~ 9 G (dotted curves). 
The long-dashed curves in the top and middle panels are 
analytic results with B ei = 3 x 10~ 7 G. In the bottom 
panel, E e g is the RS acceleration efficiency and -Ecr/^SN 
is the fraction of SN explosion energy going into CRs. The 
curve labeled 'FS' is Eqr/Esn for the FS (FS results are 
insensitive to -B e j)- In all cases, -Bism = 3 x 10~ 6 G and 
curves are truncated when p max < 2m p c. 

1000 yr. The forward shocks put ~ 50% of Esn into CRs, 
as shown by the heavy-weight solid curve labeled "FS." 6 
In the strong nonlinear regime, radical changes in the 
structure of the interaction region between the forward 
and reverse shocks occur. In Fig. 0] we compare nonlinear 
density and temperature profiles with B e] = 3 x 10~ 6 G 

6 As noted bv lBerezhko et all <l2002l) . injection may vary over 
the surface of the SNR and be significantly less where the 
magnetic field is highly oblique. They estimate that to sup- 
ply the galactic CR population the overall efficiency need only 
be ~ 20% of the maximum values obtained by DSA. 
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Fig. 4. The top panel shows the plasma density vs. radius 
and the bottom panel shows the temperature vs. radius. 
In both panels, the solid curves include the effects of ef- 
ficient DSA with a constant ejecta magnetic field B Q] = 
3 x 10~ 6 G, the dashed curves assume B ei = 3 x 1CT 8 G, 
and the dotted curves are results with no acceleration 
(i.e., test-particle). In all cases the results are calculated 
at t S NR = 1000 yr with Bi SM = 3 x 1CT 6 G. 



(solid curves) and B e j = 3 x 10~ 8 G (dashed curves) 
against profiles obtained with no acceleration (dotted 
curves). All of the profiles are plotted at £snr — 1000 yr. In 
the cases with efficient particle acceleration, the increased 
compression ratios make the interaction region narrower, 
denser, and cooler than with no acceleration. This is par- 
ticularly true for the width of the shocked ejecta region 
which, compared to the case with no acceleration, shrinks 
by a factor of ~ 5 when B C j = 3x 10~ 8 G and where the 
temperature drops from greater than 10 s K to ~ 10 6 K. 

Note that the density and temperature profiles drop 
abruptly with no apparent precursors (the small precur- 
sors seen in Fig.0J and Fig.EJbelow, are numerical from a 
finite grid). This is because our current CR- hydro model 
does not explicitly include the cosmic-ray precursor in 
the hydrodynamics. The precursor effects essential for 
the nonlinear acceleration (i.e., shock smoothing and pre- 
heating) are included in the particle acceleration model. 
For the solid curves in Fig. 0J the upstream density pre- 



cursor at the reverse shock would be ~ 0.2 pc in extent 
to be consistent with our assumption that / S k-R s k sets the 
upstream diffusion length of particles of momentum p max - 
We further note that regardless of the density precursor, 
a precursor in thermal X-rays is not expected because the 
temperature in the precursor is too low. The nonthermal 
X-ray precursor will be much narrower than the density 
precursor if the electron spectrum is limited by cooling, as 
suggested by the narrow width of the filaments. 

The changes in structure and evolution of the SNR 
brought about by efficient DSA are so large that, even 
considering the difficulties projection effects present and 
other uncertainties, they should be observable with cur- 
rent techniques. On the other hand, if radio emission is 
unambiguously observed at reverse shocks without such 
dramatic structural changes, this would be evidence for 
either magnetic field amplification beyond several /iG, or 
that DSA doesn't produce large compression ratios as the 
theory with particle escape predicts, or that the injection 
rate is considerably lower than the value ?7i n j = 10~ 3 we 
have assumed. 



4.3. Non-uniform, diluted ejecta field 

The structure of the magnetic field in SNR ejecta is clearly 
more complex than assumed in the previous section. As 
the ejecta expands, the conservation of magnetic flux will 
cause the field strength to decrease rapidly and after only 
a few years, B ei may fall below the 3 x 1CP 8 G value used 
in Fig. The highest possible initial values of the mag- 
netic field are expected for white dwarf (WD) progeni- 
tors of Typ e la supernova e and range between 10 5 and 
10 9 G (e.g.. iLiebertl Il995f) . Further enhancement of the 
WD magnetic field by convection effects could occur prior 
to explosion during the phase of quasi-static burning of 
carbon and could possibly lead to equipartition between 
the kinetic energy density and magnetic field density. If 
this occurs, fields as high as 10 10 — 10 11 G might result 
fe.g.. lRuiz-Lapuente fc Sprufflll998j) . providing the upper 
limit on the initial magnetic field intensity we consider in 
the following. 

After the explosion of the progenitor, the rapid ex- 
pansion of the ejecta will dilute the magnetic field. We 
obtain expressions for the magnetic field at time, t, in the 
core (i.e., the plateau region) of the ejecta B corc , and in 
the outer, power-law part of the ejecta Bpl, as a func- 
tion of fluid speed V , by assuming that the magnetic 
field is initially uniform in a constant density progenitor 
and then carried passively over during the explosion with 
the magnetic flux being conserved during the expansion. 
Therefore, 



V 

Bpi, = Scorn I T-r- 



(7) 

(8) 
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Here, -Bwd is the magnetic field in the white dwarf before 
the explosion, V^j is the constant fluid velocity at the core- 
power-law transition point, R a is the radius within the 
white dwarf enclosing the mass of what will become the 
core of the ejecta, and n is the ejecta density power-law 
index. Given T/^ J ax , -Esn and to ; ^ej is obtained numerically 
in the hydro simulation and 



Ra — %D 



n — 3 



3(Vniax/^ej) 



3 — n 



1/3 



0) * 



where -Rwd is the white dwarf radius before the explosion. 

Despite specifying white dwarf parameters in the 
above equations, similar behavior is expected for the ejecta 
in Type II supernovae but starting with a lower magnetic 
field at a larger radius. Basically, any pre-SN field will be 
diluted by the expanding ejecta such that B c) cx R~ 2 . 

In Fig. [SJ we show the variations of r to t, Pmax, and v c 
versus £snr for a diluted field B , as defined by eqns. 
and (El)- The frequency Vr is the critical frequency (e.g., 
iRvbicki fc Lightmanl fl9T9^ where synchrotron emission 
from an electron of p max peaks, i.e., 



v c = 3pi ax e J B/[47r(m c c) 



(10) 



Here m c is the electron mass and we have taken the sine 
of the pitch angle = 1 for convenience. 

The solid curves in Fig. [S] show the extreme case with 
-Bwd — 10 11 G. In this case, the ability of the reverse 
shock to produce ~ 10 GeV particles lasts until isNR ~ 50 
yr, after which Bo drops below ~ 10~ 8 G. The dashed 
curves show results for Bwd — 3 x 10 10 G. Despite Bwd 
being as high as 3 x 10 10 G, few GeV particles are pro- 
duced at the RS. This is contrasted by the FS (dot- 
dashed curves) where > 10 4 m p c particles are produced 
after - 100 yr. 

Even though the maximum momenta produced in 
these diluted magnetic field examples are low, large struc- 
tural changes occur early in the evolution, as shown in 
Fig. El where the density and temperature profiles are 
plotted at tsNR = 100 yr. The differences between the 
test-particle case (dotted curves) and the efficient acceler- 
ation cases (solid and dashed curves) are greater for the 
forward shocks, but are substantial at the reverse shocks. 
The structure changes at the reverse shocks lessen as isNR 
becomes greater than 100 yr. 

In Fig. we show the reverse shock synchrotron emis- 
sion predicted for the -Bwd = 3 x 10 10 and 10 11 G ex- 
amples shown in Fig. along with emission from the for- 
ward and reverse shocks where B ei is held constant at 
3 x 10~ 6 G and 3 x 10~ 8 G as labeled. For calculating 
synchrotron emission here and elsewhere, we assume the 
electron to proton density ratio at relativistic energies to 
be 0.01, similar to that observed for gal actic cosmic rays 
fsee lBaring et al.lll999HEllison et al.l l2000. for discussions 
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Fig. 5. Results from two models with diluted magnetic 
fields are shown. In all panels, the solid curves are RS 
values with £?wd = 10 11 G, the dashed curves are RS 
values with %d = 3 x 10 10 G, and the dot-dashed curves 
are FS values, which are insensitive to -Bwd- For both 
models, Bjsm = 3x 10~ 6 G and only values where p max > 
2 m p c are plotted. 



of how electrons are treated and synchrotron emission cal- 
culated in this model). These curves are calculated at 
tsNR = 1000 yr at a distance of 2 kpc, typical of SN1006. 
At 1 GHz radio frequencies, the reverse shock emission 
from even the most extreme white dwarf case with dilu- 
tion falls more than 5 orders of magnitude below that of 
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profiles with diluted ejecta fields. The dotted curves are 
results with no shock acceleration, the solid curves are 
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results with B WD = 3 x 10 10 G. 



the forward shock. 7 For the constant B c] = 3 x 10 



G 



case, the reverse shock emission is about a factor of fOO 
below the reverse shock emission with B e] = 3 x 1CP 6 G 
at 1 GHz. Note that the 3 x f (T 8 G field is too weak to 
produce X-ray synchrotro n emission as reported for SNR 
RCW 86 llRho et al.ll200l . 

The evolution of the synchrotron emission is shown in 
Fig-IHlfor our extreme diluted ejecta magnetic field £>wd = 
10 11 G. The difference in radio emission at 1 GHz between 
the forward and reverse shocks is about a factor of 20 at 
30 years and drops to more than five orders of magnitude 
at 1000 years, as shown in Fig. 



4.4. Heavy elements plasma 

In the examples we have shown so far we considered only 
acceleration in fully ionized hydrogen. However, ejecta ma- 
terial is expected to be composed mainly of heavy elements 
and, in particular, type la supernovae are essentially de- 



7 The properties of the forward shock are quite insensitive 
to the ejecta magnetic field so each of these four models have 
similar forward shock synchrotron emission. 







CD 




CO 




1 

<N 




E 






-5 


^> 




CD 




J*. 




o 




.og. 


-10 




-5 
Log 10 keV 
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void of hydrogen. There are two reasons for the accelera- 
tion process to be modified in the case of heavy elements. 

First, the acceleration time and diffusion length de- 
pend on charge. For relativistic particles, the time required 
to accelerate an ion with charge q = Ze to momentum 
p is proportional to X/Z. Likewise, the diffusion length 
of an ion with mo mentum p is proportional to 1/Z (e.g., 
tearing et alllQQflh . Therefore, for given shock parameters, 
Pmax oc Z regardless of whether p max is determined by a 
finite shock age or size. A higher p max tends to increase 
the acceleration efficiency if heavy ions are dominant com- 
pared to protons being dominant. 

Second, a species with mass number A must have mo- 
mentum > ArripC to be relativistic. Since p max only in- 
creases as Z, this tends to result in a lower fraction of 
relativistic particles, and a 7 closer to 5/3, than in the 
case where protons are dominant. This tends to decrease 
the acceleration efficiency. In a mixed plasma containing 
light and heavy ions, the modified shock structure re- 
sulting from efficient DSA produces an additional effect 
whereby ions with the largest mass to charge ratio are 
accelerated from ther mal energies most efficiently. Th is is 
described in detail in lEllison. Drurv. k, MeveJ l)l997t) but 
is not considered here where we treat only single compo- 
nent plasmas. 

In Fig. El using a constant B c] = 3 x 10~ 8 G (solid 
and dashed curves) and our extreme white dwarf diluted 
field £>wd = 10 11 G (dot-dashed and dotted curves), we 
compare r to t and p max for a hydrogen plasma and for a 
fully ionized oxygen plasma. All other parameters are the 
same as those used for the examples shown in Fig. [3] For 
the constant B e] = 3 x 10~ 8 G case, the difference in r to t 
between protons and fully ionized oxygen becomes minor 
after ~ 50 yr. For the diluted magnetic field case, however, 
there are large differences in r to t between the two species. 
The oxygen plasma shows a lower maximum r to t peaking 
at an earlier age than the proton plasma. Based on this 
alone, whatever structural and temperature changes which 
occur in the interaction region are expected to be greater 
in hydrogen dominated Type II supernovae envelopes than 
in heavy ion dominated Type la supernovae. 



5. Discussion 

5.1. Observational evidence for DSA at reverse shocks 

Projection effects and other problems make it difficult to 
reliably associate radio or X-ray nonthermal emission with 
reverse shocks in SNRs. Nevertheless, there have been re- 

ce nt claims for bo th. 

IGotthelf et al.l ( 2001^ have identified the forward and 
reverse shocks in Cassiopeia A and showed that radio and 
Si emissivity radial profiles both show a sharp rise at what 
they characterize as the reverse shock. The sharp rise and 
fairly sharp falloff in radio emissivity as one moves out- 
ward from the center of the Cas A SNR is a good indica- 
tion of the local acceleration of relativistic electrons, but, 
as noted, projection effects make the precise determination 
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Fig. 9. Reverse shock r to t's and p max 's for material con- 
taining only protons or fully ionized oxygen, as labeled. 
The ejecta magnetic field is taken to be constant at 
B e j = 3 x 10~ 8 Gin the solid and dashed curves, while a di- 
luted field with Bwd = 10 11 G is used for the dot-dashed 
and dotted curves. Only results where p m ax > 2 Am p c are 
shown. 



of the reverse shock difficult. Furthermore, IGotthelf et all 
chose to ignore a secondary Si peak near 80 arcsec 
inside the main Si peak. The main Si peak seems clearly 
associated with radio emission, but the inner secondary 
peak is not and could alternately mark the reverse shock. 

Using new VLA observations of Kepler's SNR, 
DcL anev et al ] ll2002h suggest that two distinct radio 
structures are present. These flat- and steep-spectrum 
components are partially decoupled in some areas and the 
steep-spectrum component tracks the X-ray emission seen 
by ROSAT, which is mainly line emission from shocked 
ejecta. They conclude that the flat- and steep-spectrum 
radio emission come from the forward and reverse shocks, 
respectively. However the steep regions could also mark 
the interface rather th an the reverse shock. For example, 
IVink fc Lamind |2003) interpreted the (steep) radio emis- 
sion from Cas A as arising at the interface. 
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iRho et alJ l)2002|) claim that the detailed morphology 
of soft and hard X-rays in SNR RCW 86 strongly supports 
the case for different origins for these components. They 
conclude that the hard X-rays are most likely a combina- 
tion of a synchrotron continuum plus Fe Ka emission from 
shocked ejecta material. This implies, in their estimation, 
that the reverse shock may be accelerating electrons to 
energies of the order of 50 TeV ! 

While other interpretations are certainly possible, e.g., 
relativistic electrons produced at the forward shock might 
travel to the reverse shock and brighten in radio in the 
compressed magnetic field, the above observations provide 
some evidence that reverse shocks in SNRs can directly 
produce relativistic electrons; GeV in the case of radio 
and TeV in the case of nonthermal X-rays. If this is so, 
there are a number of important consequences both for 
interpreting a wide range of astrophysical sources where 
DSA is believed to generate relativistic electrons, and for 
our understanding of the basic functioning of DSA. 

5.2. Implications 

5.2.1. Magnetic field amplification 

If relativistic electrons are accelerated at the reverse 
shock, the production of observable radio emission re- 
quires that the magnetic field upstream of the reverse 
shock be amplified by many orders of magnitude over val- 
ues expected in the expanding unshocked ejecta. If such 
field amplification is taking place , it is likely that it is 
directly associated with DSA (I.e.. Bell fe Lucekll20oi]) . 

We speculate that early in the SNR evolution (i.e., 
is nr. 10 yr) the unamplificd ejecta field may be strong 
enough to start the DSA process even with dilution. Once 
energetic particles are produced, the magnetic field may 
be amplified by them and maintained against dilution at 
a level where DSA continues. 

Shocks are widespread in astrophysics with parame- 
ters not widely different from those in SNRs. If DSA can 
amplify magnetic fields in SNR shocks, the process should 
work in other environments as well. Since the maximum 
energy individual particles obtain in DSA scales as B, 
widespread -B-field amplification will lead to a system- 
atic increase in the expected maximum proton energy pro- 
duced by astrophysical shocks. 

5.2.2. Extreme nonlinear effects 

The theory of nonlinear diffusive shock acceleration pre- 
dicts compression ratios far in excess of the test-particle 
value of rw. — 4 (e.g 



Eichlerl Il984l lEllison fe Eichlerl 



19841 : |Malkovlll998t lBelllll987t iBerezhko fc Elrisonlll999t 
BlasJl2002jr For high sonic Mach numbers, M| > M A , 



and strongly modified steady-state conditions, the com- 
pression ratio can be approximated by r to t ~ 1.3 M|^ 4 or 
by rw. ~ 1.5M 3 J 8 when Mg < M A l|Kazanas fe EllisorJ 



U.y 'tot ~ l.OiHl wneii lVlg ^ 

1986: Berezhko fe Ellison1 ll999). 



We have shown that if B Ci falls in a range 3 x 10 < 
B e j < 3 x 10~ 6 G with typical SNR parameters, r to t ~ 60 
may occur (see Fig. EJ; a value considerably larger than 
is likely to occur in any other non-radiative astrophys- 
ical environment. Compression ratios this large should 
produce unmistakable changes in the SNR structure (see 
Figs. Hand El and X -ray emission, thermal and nonther- 
mal. Confirmation of this prediction would support the 
premise that DSA is intrinsically extremely efficient. 

Unfortunately, shock compression ratios are not di- 
rectly observable for remote systems such as SNRs 
with large Mach numbers (see s ects. 2.2.3.1 and 2.4 in 
iDrurv et al.l200lilRavmondl200ll respectively). Shocks in 
the heliosphere are directly observable but have low Mach 
numbers and even here, the direct measurement of com- 
pression ratios requires multiple spacecraft simultaneously 
sampling the upstream and downstream plasmas. Until 
the launch of CLUSTER this was not possible and the 
prediction of r to t > 4 has not been clearly demonstrated in 
astrophysical shocks, alth ough indirect support does exist 
for th e Earth bow shock (Elli son. Mobius. fc Paschmannl 
Il990h . 

The large structural changes brought about in SNRs 
if ''tot 4 offer a unique opportunity to see the effects 
of extremely efficient diffusive shock acceleration. In the 
extreme cases we show here, i.e., Fig. |3 more than 90% 
of the bulk flow energy flux (in the shock rest frame) is 
placed in relativistic ions. Even in cases when r to t is not as 
extreme, acceleration efficiencies near 50% are predicted 
for both the forward and reverse shocks. Over the lifetime 
of a SNR, ~ 50% of E$n is predicted to be put into cos- 
mic rays depending on t he average injection rate over the 
surface of the SNR (e.g. lDorfilll99(A IBerezhko et"al"1ll996l 
l2002UEllison et alJl20o3) . The energy which goes into rel- 
ativistic ions comes out of the bulk thermal plasma and 
produces a drastic reduction in the shock temperature. 



5.2.3. Cosmic-ray production 

Direct observations of shocks in the heliosphere and most 
theories of DSA show that collisionlcss shocks put far more 
energy into ions than electrons. Thus, even though the 
presence of relativistic electrons suggested our descrip- 
tion of strong nonlinear effects, the signature of nonlin- 
ear DSA in the structure and evolution of the radio and 
X-ray emitting interaction region between the forward 
and reverse shocks, will be evidence for the efficient pro- 
duction of cosmic-ray ions, not neces sarily electrons (see 
IBerezhko. Ksenofontov. fc VolldEooll for a discussion of 
SN 1006 in this regard). Detection of pion-decay 7-rays 
would be a more direct confirmation that relativistic ions 
are produced in SNRs, however this has not yet been un- 
ambiguously done. It is also clear from 7-ray models of 
young SNRs, that parameters can be chosen where effi- 
cient DSA occurs, but either the 7-ray flux is below de- 
tectable levels, or detectable TeV 7-rays are dominated 
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by inverse- C ompton emission from electrons rather than 
protons (e.g.. lEllison. Slane. fc Gaenslerll20o"l|) . 

Modeling the SNR structure offers two advantages. 
First, if r to t > 4 is inferred from observations, it is evi- 
dence for the efficient production of CR ions whether or 
not a pion-decay feature is observed since the observability 
of pion-decay 7-ray's depends on other factors besides ac- 
celeration efficiency (e.g., ambient density). Second, if the 
structure is inconsistent with efficient DSA and r to t & 4, 
this is clear evidence that the efficient shock acceleration 
of ions is not occurring. Due to the freedom of 7-ray mod- 
els, the lack of a 7-ray detection is unlikely ever to be able 
to eliminate the possibility that efficient DSA is occurring. 

5.2.4. X-ray emission 

For efficient particle acceleration, the postshock densities 
are larger and the postshock temperatures smaller than 
in the test-particle case (see Figs. 0] and EJ). As a con- 
sequence, the heating of electrons in the downstream re- 
gion by Coulomb interaction with the population of pro- 
tons (T p ~ 1836 T e without further heating at the shock) 
will be more efficient than in the test-particle case. In 
the shocked ejecta, the electron temperature can be al- 
most equal to that of the protons for an injection of 
~ 10~ 2 , while in the shocked ambient medium the elec- 
trons may reach up to 30% of the proton temperature 
l|Decourchelle fe Ellisonll200l|) . 

Another constraint on the efficiency of particle acceler- 
ation comes from the observation of a strong Fe K-alpha 
line in the shocked ejecta of young SNRs (like Cas A, 
Kepler and Tycho). In Kepler, for efficient particle accel- 
eration at the reverse shock, it was shown that the shocked 
ejecta temperature gets t oo low to produce the Fe K-line 
l|Decourchelle et all200o|) . Line excitation by the nonther- 
mal population can be invoked, however the ionization 
state of iron is expected to be very low even when tak- 
ing into accoun t the ionization from a nonthermal power- 
law population i Porauet et all200~iL Thus, the prediction 
from efficient DSA that postshock temperatures are low 
in the RS presents a problem in SNRs where a strong Fe 
K-alpha line is observed. 



5.2.5. Ionization fraction of ejecta material 

For DSA acceleration to occur at all, the unshocked ejecta 
material must be fully ionized or, at least, have a siz- 
able ionization fraction. Otherwise the magnetic turbu- 
lence necessary to scatter particles will be damped (e.g., 
iDrurv et al.l ll9961. The expanding ejecta will cool rapidly 
and would be largely neutral unless ionized by some 
source. This could possibly be X-ray emission from shock 
heated gas or the cosmic rays may contribute to the ioniza- 
tion themselves. If the ionization fraction is initially large 
enough for some cosmic-ray production to occur, these 
CRs may further ionize the precursor material, making 
acceleration more efficient, etc. 



The Balmer-dominated s pectra of nonradiative sh ocks 
in a number of SNRs (e.g.. iGhavamian et al.ll200ll) has 
been interpreted as charge exchange between protons and 
neutrals. If synchrotron radio and/or X-ray emission is 
observed at these shocks as well, which is possibly the case 
in the northeastern rim of Tycho, this would indicate that 
DSA is at work even in a partially neutral medium. The 
reported presence of reverse shock radio emission suggests 
that the ejecta material is ionized, at least in the region 
upstream of the shock, where DSA is taking place. 

6. Conclusions 

It has long been believed that forward shocks in SNRs 
sweep up and accelerate ISM ions to produce galactic cos- 
mic rays and to accelerate electrons to produce in situ ra- 
dio emission. The mechanism most likely responsible for 
this is DSA, which is predicted to be extremely efficient, 
i.e., > 50% of the ram kinetic energy may go into rela- 
tivistic ions. At reverse shocks in young SNRs, however, 
straightforward estimates of the diluted ejecta magnetic 
field from the pre-SN white dwarf or massive star show 
it to be many orders of magnitude below that required 
to either accelerate relativistic electrons from the ther- 
mal background, or to produce observable radio intensities 
from background cosmic-ray electrons. 

Nevertheless, some recent observations provide evi- 
dence for radio emission associated with reverse shocks 
in two SNRs: Kepler and Cas A, and X-ray synchrotron 
continuum emission in RCW 86. If the reverse shocks in 
these remnants are accelerating electrons to GeV or even 
TeV energies via DSA, it immediately suggests that the 
magnetic field at the reverse shock is orders of magnitude 
higher than expected to produce radio and higher still to 
produce X-ray synchrotron. If true, the most likely expla- 
nation is that the acceler ation process is amplifying the 
magnetic field, perhaps as lBell fc Lucekl l|200ll) have sug- 
gested. The importance of shock acceleration in a wide va- 
riety of astrophysical objects, and the strong dependence 
of DSA on the magnetic field, make it critically important 
to first, verify the difficult reverse shock observations and 
second, to explore the ramifications of efficient DSA at 
reverse shocks in young SNRs. 

While radio observations imply magnetic field ampli- 
fication to values far larger than the diluted progenitor 
field, there is, as yet, no way to precisely determine the 
value. This opens the possibility that the ejecta field B c] 
is large enough to produce observable radio emission, but 
still far lower than normal ISM values. We have shown 
here that, for typical SNR parameters, having 3 x 10 -8 < 
B e j < 3 x 10~ 6 G results in extremely large nonlinear ef- 
fects in DSA, i.e., compression ratios ^ 4 and shocked 
temperatures <C than test-particle temperatures. If these 
nonlinear effects occur, they will produce large changes in 
the structure and evolution of SNRs which should be ob- 
servable with existing instruments. We have detailed these 
effects in the remnant hydrodynamics and estimated syn- 
chrotron emission for a limited range of parameters. 
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Supernova remnants may offer the best known labo- 
ratory for studying both magnetic-field amplification and 
DSA. Current and future ground and space-based obser- 
vatories offer high spatial and energy resolution of sev- 
eral remnants and provide information that is available 
nowhere else. Shocks in the Heliosphere are accessible to 
spacecraft and a great deal has been learned of their prop- 
erties. However, the point-like nature of heliospheric ob- 
servations, the low Mach number and small size of helio- 
spheric shocks which limits particle acceleration, and the 
unique geometry of the intensely studied Earth bow shock, 
limit what can be learned and transferred to other astro- 
physical systems. The difficult plasma physics has also lim- 
ited the success of analytic investigations, and direct PIC 
computer simulations are decades away from being able 
to simulate the injection and acceleration of a electron- 
proton plasma to relativistic energies. 
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